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Abstract 

We study implications of recent data on neutrino mixing from T2K, MINOS, Double Chooz and fi — )• 
67 from MEG for the Zee model. The simplest version of this model has been shown to be ruled out 
by experimental data some time ago. The general Zee model is still consistent with recent data. We 
demonstrate this with a constrained Zee model based on naturalness consideration. In this constrained 
model, only inverted mass hierarchy for neutrino masses is allowed, and ^13 must be non-zero in order to 
have correct ratio for neutrino mass-squared differences and for mixing in solar and atmospherical neutrino 
oscillations. The best-fit value of our model for ^13 is 8.91° from T2K and MINOS data, very close to 
the central value obtained by Double Chooz experiment. There are solutions with non-zero CP violation 
with the Jarlskog parameter predicted in the range ±0.039, ±0.044 and ±0.048 respectively for a la, 2a 
and 3a ranges of other input parameters. However, without any constraint on the 0i3-parameter above 
respective ranges become ±0.049, ±0.053 and ±0.056. We analyse different cases to obtain a branching 
ratio for — )• 67 close to the recent MEG bound. We also discuss other radiative as well as the charged 
trilepton flavour violating decay modes of the r-lepton. 
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I. INTRODUCTION 



The standard model (SM) of electroweak and strong interactions based on the gauge group 
SU(3)c X SU{2)i X U{1)y has been rigorously tested in many ways, in different sectors involving 
gauge bosons, quarks and leptons. The flavour physics in the quark sector including CP violation 
is well described by the Cabbibo-Kobayashi-Maskawa (CKM) mixing matrix in the charged current 
interactions of the W boson with three generations of quarks. Corresponding interactions in the 
leptonic sector are not on the same footing. In the simplest version of the SM, there are no right 
handed neutrinos and the left handed neutrinos are massless. This theory conserves lepton number 
of each family separately, and no flavour changing neutral current (FCNC) processes are allowed in 
the leptonic sector. However, the observation of neutrino oscillations requires both neutrino masses 
and flavour changing processes. 

Study of neutrino oscillations in the past few years has led to an insight regarding mass differences 
and mixing between the three light neutrinos. The SM can be easily extended to give neutrino 
masses and explain their mixing through a mixing matrix in the charged lepton currents analogous 
to the CKM matrix in the quark sector as was first suggested by Pontecorvo-Maki-Nakagawa- 
Sakata (PMNS) [1]. Experimentally, not all the parameters of this matrix are known; in particular 
whether CP is violated in the leptonic sector. Further understanding of the smallness of neutrino 
masses requires new physics, such as loop induced masses as in the Zee model [2], or a seesaw 
mechanism which can be implemented in several ways, referred to as Type I [3] , II [4] , III [5] seesaw 
models. These models, besides generating neutrino masses and their mixing, also induce other 
FCNC processes in the leptonic sector causing new lepton flavour violation (LFV) phenomena. 
There are strong constraints on LFV interaction [6]. The Zee model has been extensively studied 
in the literature [7]. Data from T2K [8] and MINOS [9] a few months ago provided some new 
information on the mixing angles in Vpmns that the last mixing angle ^13 is non-zero at more 
than 2a level. The combined data analysis gives the confidence level at more than 3cr which put 
the well known tri-bimaximal mixing pattern into question. Very recent data from Double Chooz 
[10] also indicate that may be non-zero. Models for neutrino masses and mixing are therefore 
further constrained. The hint of non-zero ^13 sparks a few analysis in various ways, see for example 
[11]. In most of the models for neutrino masses and mixing, there are LFV interactions. For a 
study of testing a model, it is therefore also necessary to consider LFV processes. Recently MEG 
collaboration [12] has report a new upper bound of 2.4 x 10~^^ for yU — 67 branching ratio which is 
about 5 times better than previous one [6]. This improved upper bound may have implications on 
models for neutrino mass and mixing. In this work we confront the Zee model with the new data 
on neutrino mixing and different LFV processes. 

A. The Zee model 

Neutrinos being electrically neutral allows the possibility that they are Majorana particles. There 
are many ways to realize Majorana neutrino masses. Even with the restrictive condition of renor- 
malizability, there are different type of models. Without introducing right-handed neutrino N^, 
one can generate Majorana neutrino masses by introducing Higgs representations. With different 
Higgs representations, one can generate neutrino masses at the tree or loop levels. The Zee model 
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is a very economic model for loop induced neutrino masses providing some reasons why neutrino 
masses are so much smaller than their charged lepton partners. 

In the Zee model, in addition to the minimal SM without right-handed neutrinos, there is another 
Higgs scalar doublet representation $2 : (1? 2, 1) beside the doublet $ = $1 already in the minimal 
SM, and a charged scalar singlet /i^ : (1, 1,2). Terms in the Lagrangian relevant to lepton masses 
are 



L = + fl<^2)ER - Llfia^LLh^ - ii^lia^^h^ + h.c. (1) 

where = CL^ is the charge conjugated : (1,2,-1), Er : (1,1,-2) is the 5'f/(2)— singlet 
right-handed charged lepton and /i is the tri-scalar coupling constant with mass-dimension one. 
The Yukawa coupling matrices fi are arbitrary while / is anti-symmetric in exchanging the family 
indices. 

This leads to a mass matrix for the charged leptons 

mE = ^ (sin/3 f| + cos/3f|) , (2) 
V 2 

where v = ^^/vf + V2 and tan/3 = vi/v2 with 

We will work in the basis where charged-lepton mass matrix, m^, is diagonalised, such that = 
diag (me, m^,m^). 

The would-be Goldstone bosons and z "eaten" by the and Z bosons and the other two 
physical components, and a are given by 



_ f cos/3 — sin/3\ / f^\_f cos/3 —sin/3 



w+ J I sin/3 cos/3 J J ' v) W^l^ cos/3 ) U)% 



(4) 



For the two CP even fields 0_ri,_r2) we will use their mass eigen-states H and h linear combinations, 



H\ _ f cosa — sina\ / </)5ji 
h j I sina cosa j \4>\2 



(5) 



There will be mixing between and if is not zero. Without loss of generality, one can 
write the mass eigen-states and as 



_ f cos^z — sin^z\ / h 
h+J ^ ysine, cos9, J \H+ J ' 



(6) 



We have the following lepton-scalar couplings 
L = -ul 
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-2z/£/EL(cos^,h+ + sin^^h+) 
^ vsin/3 

I vsin/3 

yvtanfj V2sin/3 

The neutrino mass matrix, defined by (l/2)i/£Mj,z/2,, is related to the model parameters as 



ZTo — sin(p — a 

- ( mucosa 
-Ll —5 cos(p — a) 




(7) 



(8) 



where t/ = Vpmats is the PMNS mixing matrix, and 

A = sin(2^,) In (ml+/ml+^ / 8V2iT\tan(3 . (9) 

In the simplest Zee model, /2 = 0. It can be easily seen from the expression for M,^ that in this 
case, the resulting mass matrix has all diagonal entries to be zero. This type of mass matrix has 
been shown to be ruled out by data [13]. This is basically because that it cannot simultaneously 
have solution for \Vf^s\ close to l/\/2 and \Ve2\ close to as data require. With /2 non-zero 

it can fit experimental data. In this case, however, we encounter a problem which is common to 
many models beyond the SM: that there are too many new parameters. Additional theoretical 
considerations help to narrow down the parameter space. To this end, it has been proposed that an 
interesting mass matrix can result if one imposes the requirement that no large hierarchies among 
the new couplings, that is, all f^^ and f2 are of the same order of magnitude, respectively, from 
naturalness consideration. From the expression for the neutrino mass matrix one sees that all terms 
are either proportional to charge lepton mass m/ or mass-squared mf. Since >> m^,me, the 
leading contributions to the neutrino mass matrix are proportional to f^'^m'^ and /J^m^. To this 
order, one can write the neutrino mass matrix using the five input parameters x, j/, z, a and 5 only, 
without loss of generality, as 

/ 1 {ye'^ + x)/2 z \ 

M^ = ai {ye'^ + x)/2 xye'^ xz . (10) 
y 2; xz J 

Here a is the absolute value of the 11 entry Mn, 

Mn = -2Avm^rf^ycos^, (11) 

of and z is the absolute value of 13 entry M13 divided by a with 

Mrs = Arm,{mr-vjY/cosf3), (12) 

^ = I /'^'^ I / 1 Z*^"^ I is the absolute value of the ratio of M23 to M13 and y = \M22\/xa with 

M22 = -2Avmrf''^fl^/cosP. (13) 

One can choose a convention where all the above parameters are real except /J'^. 
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II. NEUTRINO MASSES AND MIXING 



The mass matrix M^, is of rank two with the two eigen masses given by 



— (|1 + xye^'Y + (1 + x'){l + y' + Az') 
±2|1 + xye^^l V(l + + y'^ + 4z^)^ . 



(14) 



If one identifies mi 
given by 



m_, 7712 = 1^+ and = 0, the corresponding Vpmns niixing matrix is 



V, 



PMNS 



V2Ni 
V2Ni 
Ni 



V2N2 
xB+(l+xye'^)ye-"^ 

V2N2 
V2{l+xye^^)z 
N2 



2X2 \ 

2z 

ye^^-x 
N3 / 



(15) 



where 



B 



xye'^ I ^/{l + y^ + 4:Z^)/{l + x^) 



\l + xye'^\ (^(1 +y^ + Az^) t |1 + xye'^\^y{l + y^ + Az^)/{1 + x^) 
{{l + x^){l + y' + 4z^)-\l + xye''\Y' . 



1/2 



(16) 



Other identification would imply exchanges of column vectors in the above mixing matrix. 
The non-zero masses in the above basis have Majorana phases. They are 



7n± 



-(1 



xye 



(^1 ± + a;2)(l + y'^ + 4z^)/\l + xye^^\'^^ 



(17) 



Since ma = 0, the Majorana phase matrix P can be written as P = diag(l,e*°, 1). The simple 
structure of our model allows us to determine the Majorana phase a analytically. We defined the Ma- 
jorana phase with the standard Particle Data Group (PDG) form for mixing matrix Vpmns [6], where 
the el- and e2-elements are real. With the help of the equation, = (1^mns-P)*"^*^'^^(^mns-P)^ 
one can have the following relation 



2a 



With the fact that the mass ratio m2/mi is real and negative, we obtain 

2xysmS{l + x^)(l + xycos(5) 



cot a 



(y2 — + 4^2) (1 + x'^ + y'^ + 2xycos6) 



(19) 



The Jarlskog parameter J = Im{yeiV*2V*iVfj_2) is given by 

2xyz'^ y/l + x^smS 



J 



((1 + x2)(l + y2 + 4z2) - |1 + xyei'5|2)|l + xyei^|v/(l + y2 + 4z2) 
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2(ml 



(20) 



This model only allows inverted neutrino mass hierarchy solution with mi and m2 identified 
with m_ and m+ in the above, and ms = 0. For normal neutrino mass hierarchy, that is mi = 0, 
one would require the 13 entry, xz/N^ = iKil, to be of order a/2/3 and the 12 entry, |(i? + (1 + 
xye^^) / \/2N2\ = iKsl, to be small in eq.(15). Combining with the information that Am'^i/ = 
m^/(m^ — m^) is small, we find that there is no solution. 

With known ranges for |V^2|, l^'^sl? and Am2i/Am|2, even without knowing the value for Ves, it 
is constrained to be non-zero at more than la level. This model requires a non-zero Vgs in consistent 
with T2K, MINOS and Double Chooz data. 



Parameter 


W/lO-^eV^ 


AmVlQ-^eV^ 




sin26'23 




Best-fit 


7.58 


2.35 


0.312 


0.42 


0.025 


la range 


7.32 - 7.80 


2.26 - 2.47 


0.296 - 0.329 


0.39 - 0.50 


0.018 - 0.032 


2(7 range 


7.16 - 7.99 


2.17 - 2.57 


0.280 - 0.347 


0.36 - 0.60 


0.012 - 0.041 


3a range 


6.99 - 8.18 


2.06 - 2.67 


0.265 - 0.364 


0.34 - 0.64 


0.005 - 0.050 



TABLE I: Ranges for mixing parameters obtained in Ref.[14]. The two mass-square differences are defined 
as 5m? = 7712 — ml and Am^ = mg — (m^ -|- m^) /2 



Evidence of non-zero reactor angle is published by T2K. When these data are combined with the 
data from MINOS and other experiments it clearly indicates a large deviation of the reactor angle 
form zero value [14] [15]. In the following we present our results for allowed parameter space using the 
combined neutrino mixing data, including the recent T2K and MINOS results, given in Table-I from 
Ref.[14] with the new reactor fiux estimate. The new Double Chooz data, sin^(2^i3) = 0.085 ±0.051 
is consistent with the ranges given in Table-I. 

The best-fit values of the mass matrix parameters are: 

X = 0.255, y = 4.100, z = 1.790, a = 0.017 eV, S = 180°, (21) 

and, the corresponding output for the mixing angles and mass-squared differences are given as 
follows 

sin^^is = 0.3163, sin2^23 = 0.4033, sin^^ia = 0.0256, 

6m^ = 7.51 X 10"^ eV^ Am^ = -2.36 x 10~^ eVl (22) 

The fact that our solutions are for inverted neutrino mass hierarchy case only is indicated by the 
negative sign on Am2 (I^I±!I^). In this case, solutions are in well agreement with la range 

of experimental data of Table-I. However, there is no solution to be consistent with the normal 
hierarchy case. Using the T2K and MINOS data above the best-fit value of our model for ^13 is 
8.91°, very close to the central value obtained by Double Chooz experiment. 

Variation of input parameters x, y, z, a and S that satisfy the experimental neutrino data is shown 
in Fig.-l. In first row x vs y, x vs z and y vs z are respectively shown in first, second and third plots. 
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Input parameters are bounded to be within the larger-grey (purple-circle) area to generate neutrino 
data within 3cr ranges. To satisfy neutrino data within 2a and la ranges, the input parameters 
have to be within the mid-white (yellow-triangle) and smaller-dark (red-star) areas respectively. 
We see here that the allowed range of the x-parameter lies in between (0.23, 0.29), (0.19, 0.34) and 
(0.15, 0.40) for lo", 2a and 3(T respectively. Corresponding limit on y-parameter one can read as 
(4.0, 6.4), (3.4, 13.6), (3.2, 20.0) and that for the z-parameter is (1.9, 2.8), (1.6, 5.2), (1.0, 8.0) 
respectively. 



iiiiiiiiiiiiiiiiiiiiiiiiii 








14 0.21 0,28 0.35 0,42 0.14 0.21 0.28 0.35 0.42 4 8 12 16 20 24 




«5 1 00 ~- 



FIG. 1: Variation of different input parameters (x, y, z, a, (5 ) are shown here. Input parameters are bounded 
to be within the larger-grey (purple-circle) area to generate neutrino data within Za range. To have neutrino 
data within 2a and Icr ranges, the input parameters will have to be within the mid-white (yellow-triangle) and 
smaller-dark (red-star) areas respectively. 

On the second and third row, respectively, we have shown the variation of the parame- 
ters a and 5 with the x, y and ^-parameters. Here, we see that the parameter a can be 
in between (0.015 eV, 0.019 eV), (0.011 eV, 0.020 eV) and (0.007 eV, 0.023 eV) for la, 
2a and 3a respectively. The same for the CP-phase parameter 5 is respectively given by 
(140°, 210°), (110°, 230°) and (100° - 250°) for la, 2a and 3a with a central value at 180°. 
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The allowed range for x (= |/'^^|/|/'^^|) is a good evidence that our solutions are consistent with 
our naturalness assumption that the non-zero fij should be the same order of magnitude. 

The (1,1) entry m^^^ of the neutrino mass matrix My can induce neutrinoless double beta decay. 
In our model m^^^ = a is not zero, neutrinoless double beta decay can, therefore, happen. Experi- 
mentally, nip^^ is constrained to be ^2 eV [6]. In our analysis we see in Fig.-l that the parameter 
a is allowed upto 0.019 eV, 0.020 eV and 0.023 eV for la, 2a and Scr cases respectively. Thus, our 
model is well below the allowed range for neutrinoless double beta decay effective neutrino mass 
parameter. 




FIG. 2: Variation of the Jarlskog parameter with different input parameters (x, y, z, 5 ) are shown here. In 
this case neutrino mass-square differences as well as all three mixing angles, including ^13, are within the ranges 
set as given in Table-I. 

There are some ongoing and planned experiments, for example GERDA [16] and Majorana [17], 
to search for the neutrinoless double beta decay evidence. Both of these experiments are mainly 
intended to lower down the upper limit, < 2.2 eV, on the effective mass parameter in the neutrinoless 
double beta experiment upto 0.1 — 0.3 eV. These limits are much higher than the allowed range of 
the parameter "a", ~ 0.02 eV, in our analysis. Thus, it will remain inconclusive in terms of the 
upper limit. 
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However, there is also a lower limit on the neutrinoless double beta decay effective mass pa- 
rameter from the Heidelberg-Moscow experiment [18], 0.1 eV^ rriee ^ 0.56 eV at (95% C.L.). 
Once one include a ± 50% uncertainty of the nuclear matrix elements, the above ranges widen 
to 0.05 eV^ iriee ^ 0.84 eV at (95% C.L.). The lower limit is much higher than the predicted range 
of "a" from our analysis. The main concern is that this result, the lower limit, is not confirmed by 
any other independent experiment. Here, GERDA in very near future during it's phase-I will test 
the result claimed by the Heidelberg- Moscow experiment. This will be crucial for our discussion in 
terms of it's lower limit. Our model parameters have to be reconsidered once the lower limit set by 
the Heidelberg-Moscow experiment is confirmed by GERDA during it's phase-I run. 

Consequence of a nonzero input phase 6 will appear in the output Dirac CP-phase parameter. 
This is translated into a significant deviation of the Jarlskog parameter from zero value. We have 
shown the variation of the Jarlskog parameter in Fig.-2 with respect to the input parameters x, y, 
z and S. We see that the allowed ranges of the Jarlskog parameter in our scenario lies in between 
±0.039, ±0.044 and ±0.048 respectively for a la, 2a and 3a. 




FIG. 3: Variation of the output (left) and input (right) CP-phase with respect to the sin^^ia. Here, sin^^is 
is set to be free parameter while all other neutrino mixing angles as well as the mass-square differences are 
constrained. 
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One important outcome of our discussion is that the reactor angle is constrained to be nonzero. 
To have a clear picture, we have shown the variation of the sin^^is with respect to both the input 
and output phase angle in Fig.-3. In this figure sin^^is is free from any constraint while all other 
neutrino mixing angles as well as the mass-square differences are constrained to be within la, 2ct 
and 3a ranges as given in Table-I. The variation of the output (input) CP-phase angle with sin^^is 
is shown on the left (right) panel in the figure. Here, we see that the lower limit on the reactor angle 
6*13 is clearly separated to be nonzero. The corresponding lower limit on the sin^^is are 0.024, 0.016 
and 0.011 for la, 2a and 3a cases. The best-fit value predicted from our analysis for ^13 comes out 
as 8.91° which is very close to the central value obtained by Double Chooz experiment. In the same 
analysis, the Jarlskog parameter is constrained to be in between ±0.049, ±0.053 and ±0.056. 

Finally, we consider the variation of the Majorana phase (a), defined in eqn.(19), in terms of 
the input parameters x, y, z and 6. As before, here, we consider the variation in such a way that 
the input parameters lie within a range that will generate neutrino data in la, 2a and 3a as stated 
earlier. These three different respective zones are shown in Fig.-4, with same notation used in 
previous figures. The central value of the Majorana phase is around —90°. 




0.14 0.21 0.28 0.35 0.424 6 8 10 12 14 16 18 20 22 24 

X y 




I I I I I I I I I I I I I I I I I I I I I I I 1^1 I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



10 120 150 180 210 240 270 
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-90 
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-92 
;93 



FIG. 4: Variation of the Majorana phase with different input parameters (x, y, z, 5 ) are shown here. The 
input parameters are within la, 2a and 3a ranges (see text). 
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III. RADIATIVE FLAVOUR VIOLATING DECAYS ^ 67 AND r ^ 6(^)7) 



In previous section, we have obtained allowed ranges for different input parameters that satisfy 
the current neutrino data. It is desirable to check whether these parameter space are compatible 
with other phenomenological aspects. From eqn.(7), we see that non-zero elements of two Yukawa 
coupling matrices /2 and / contribute to the LFV interactions mediated, respectively, via charged 
{hi^h^) and neutral {H,h,a) scalars as shown in Fig.-5. In this section we discuss a few such 
radiative decay modes, like, the yU — > 67 or r — )■ fj^{e)'~f. Absence of any such lepton flavour violating 
process in SM, this analysis can impose stringent constraints on the parameters space. 



/ \ 



/ \ 
/ \ 
1 . , I 



FIG. 5: Diagrams contributing to the /x — > 67 process at the one loop level. 



The matrix element M for the decay mode /j — > /j7 is given by, 

M{U ^ /,7) = ^efjiCi'PL + C'^PR)a,,helq, , (23) 

here, e, being the photon polarisation vector and a^^ = z{7^,7^}. It corresponds to an effective 
Lagrangian L^ff of the form 

Leff = IKClPl + CRPR)a^^fF,, , (24) 

where, we have suppressed the flavour indices / = (e, /i, rY and C matrices are defined as Cl{r) = 
(^L{R)) with Cr = while F^y is the field strength tensor of the photon field. 

One can easily read /i — 67 decay width [19] as 

r(/i^e7) = '^{\C,\' + \Cr?). (25) 

and, in our model one can have the explicit expression for Cl and Cr in Appendix-A. 

Like other areas of new physics, here as well, large number of free parameters is a drawback 
of the model. So, to have a better understanding on how various parameters affect the branching 
ratios, we will categorically discuss them as needed. Certainly, our interest is to find that range of 
parameter space which can lead to the new experimental upper bound, Br(/i — e7) < 2.4 x 10~^^, 
set by the MEG collaboration [12] . 

In the following analysis we will set the input parameters x, y, z, a and 6 at their best-fit values 
given eqn.(21), for illustrations. Once we fix these parameters at their best-fit values, all the 
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elements of eqn.(lO) have definite values. We can determine different elements of the / and /2 
couplings in terms of the parameter A, defined in eqn.(9), using the relations from eqn. (11-13) once 
we choose /J"^ equals to zero. So, in the following discussion, to solve for other coupling constants, 
we will always set /J"^ = 0, and also all other elements of f^^ and f2 not appearing in My set to he 
to zero. The relevant nonzero elements of the coupling constant matrix / are Z^"^, Z'^'^, hence due to 
anti-symmetry /'^'^, /'^'^, while the nonzero elements of /2 are /J"^, /J'^ only. 

In this case, hf contribution to /i — )■ 67 only comes from non-zero f, while for r Ai(e)7, the 
contribution comes from non-zero /2 only. On the otherhand, three neutral scalars i/, h and a 
contribute both to /i — 67 and r — /u(e)7 from non-zero /2. We discuss our results for different 
cases when either of the coupling constants /2 or / has a dominant contribution and when both the 
couplings / and /2 have comparable contribution to the branching ratio Br{^ — )■ 67). 

1. When either f2 or f dominate the contribution 

In this section, we will discuss how the branching ratio will depend on the Yukawa couplings 
when either, /2, associated to the neutral Higgs scalars or, /, associated to the charged Higgs scalars 
has significant contribution. Nonzero elements are derived in terms of the parameter A. However, 
this parameter itself and hence the branching ratio for fi ^ e'-f depends on different mixing angles 
as well as on the mass-squared ratio of two charged scalars. Here, for illustrations, we have chosen a 
set of values of the mixing angles a = n/6 and 6^ = vr/lO while we can read tan/3 as the parameter 
A varies, once we fix the charged scalar mass ratio. We have considered the neutral Higgs scalars 
masses to be free parameters and their values chosen here are within the experimental limit. As 
an example, we set the neutral scalars h, H and a masses at 200 GeV, 250 GeV and 300 GeV 
respectively. 

In Fig.-6 we have shown the variation of the parameter A for different values of (M^+/M^+) 
starting from 1.2 to 2.0 with = 250 GeV. Once, the parameter A is chosen, the coupling con- 
stant /^"^ is determined through eqn. (12), and hence, f^"^ through the input value of x— parameter. 
This will fix the parameters /™ and /g'^ through eqn. (11) and eqn. (13) respectively. Dependency 
on different mixing angles is discussed in Fig.-7. 



In Fig.-6, to have the desired branching ratio Br{fi 67) ~ 10^^^, we see that there are two 
clearly distinct regions of the parameter A - one for a lower scale of A ~ 10~^GeV~^ and other in 
a much higher scale, A ~ lO'^^GeV"^. For the left hand region that means for a relatively lower 
region of A, the coupling constants / are much larger in comparison to /2. For example, with the 
mass ratio equals to 1.4 ( red-dash line) we have the branching ratio Br{fi — )■ 67) = 2.39 x 10~^^ 
for a value of A = 3.87 x lO^^'^GeV"^. With this value of A, the / coupling constants are given 
by = 0.0249 and /^^ = 6.34 x 10"'^. However, corresponding /2 couplings are much smaller 
and given by /J^ = —5.43 x 10"® and f2'^ = (1.33 -|- i 1.78) x 10~^. We see in this lower range, 
the branching ratio decreases with the increase of A what we can intuitively understand as follows. 
Note that, here in the lower range of A, the main contribution comes from the charged Higgs-sector 
only, as / >> /2. From eqn. (11-13), it is evident that for a fixed set of input parameters the product 
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FIG. 6: Variations of BR(^ — )■ 67) with the parameter "Log A" for different values {M^+ /Mf^+) ranging from 



1.2 to 2.0 with 



250 GeV 



of / and A is constant. So, an increment of the parameter A reduce the couphng constant / and 
hence the branching ratio. 

On the other hand, for a comparatively large value of A, ~ 10~^GeV^^, the coupling constants 
/2 has much larger contribution to that of /. For example, with the charged scalar mass ratio 
equals to 1.4 ( red-dash line) we have the branching ratio -Br(/i — 67) = 2.21 x 10~^^ for a value 
oi A = 6.43 X lO-^GeV"^ In this case, corresponding different /2 coupling constants are given by 
fl<^ = -1.97 X 10"^ and f^^ = (4.83 + i 6.47) x 10-^ while, the much smaller values for the / 
couplings are given as /^"^ = 1.50 x 10"'^, /^^ = 3.83 x 10"^. 

In this region, the characteristic of the graph is opposite to that at the lower range. The 
main contribution to the branching ratio comes from /2 only. From eqn. (11-13), we see that /2 is 
independent of the parameter A. One may, thus, expect a constant value of the branching ratio. 
However, from eqn. (9), we see the tan/3 is inversely proportional to A. Hence, an increment of A 
implies the increase of the factor l/tan/3, which is appearing in the graphs. 



Since the parameter A depends on the mixing angles as well, to see their impact on the branching 
ratio, considering the charged scalar mass ratio equals to 1.5, we have shown the variation of 
BR(/i — > 67) with tan/3 and the mixing angle 6^ in Fig.-7. We see that in either case the branching 
ratios saturate at their lowest values as the arguments increase. 

It is interesting to check if the set of parameters obtained for the /i — )■ 67 channel will satisfy other 
LFV radiative decay processes, namely, r — )■ /i7 or r — )■ 67 processes. The current experimental limit 
on the branching ratios of these two decay channels are i?r(r — )■ fi'-f) < 2.5 x 10"''^ and Br{T — )■ 67) < 
1.8 X 10-^ [6]. 

In the lower region of A in Fig.-6, where f'^'^, /'^'^, /'^'^ and f'^'^ only significant ones, while /2 is 
negligible, there will not be any significant contribution to the r — )■ n{e)'-f branching ratios. However, 
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FIG. 7: Variations of BR(/i — 67) with different mixing angles in the charged and neutral Higgs sectors 

had we experimentally observed any of these decay modes we must then analyse the model with 
corresponding non-zero coupling constants. 

On the other hand, for a larger range of A, ~ 10~"^GeV~^, the coupling /2 only dominantly 
contributes to the branching ratio. Here, /J^ and /2 will open a channel for r to radiatively decay 
into /i(e)7. In Fig.-8, we have shown the variation of the branching ratios with the parameter A for 
the larger range only. Considering the same set of input parameters we see that the corresponding 
r-decay branching ratios can reach to the experimental limit with increase of A for different values 
of charged scalar mass-ratios. 




2e-05 4e-05 6e-05 8e-05 0.0001 0.00012 0.a)OI4 0.00016 0.00018 
A (GeV ') 




2e-05 4e-05 6e-05 8e-05 0.0001 0.(K)012 0.00014 0.00016 0.00018 
A(G<!V') 



FIG. 8: Variations of BR(r — fi^) (left) and BR(t — )■ 67) with the parameter "A" for different set of values 
(M^+/M^+) ranging from 1.2 to 2.0. 
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2. When both f and f2 significantly contribute 



Here, we study the parameter space when both the couphngs / and /2 will have significant 
contribution to the branching ratio. The parameter space chosen in previous discussion can not 
lead to MEG result. However, we realise that almost degenerate charged Higgs masses can fuUfil 
the criterion. Variation of the corresponding branching ratio BR{fi — 67) is shown on the left 
panel of Fig.-9 for different values of the charged Higgs scalar mass ratio (M^+/M^+) ranging from 
1.00002 to 1.00010. In this case, other mixing angles as well as the neutral Higgs masses are kept 
unchanged as like before. 

To understand the characteristic of these graphs, we note that in higher range of A, the couplings 
/2 is large while / is negligible, and for A of the order of 10~^°GeV~^ the case is opposite. It is 
obvious that for the range of A ~ (9(10~^)GeV~^, which is almost half-way from two of the above 
mentioned scale of A in Sec.-III.l, would give same order of values for both the couplings / and 
/2. However, corresponding branching raito is much smaller in magnitude compared to the MEG 
data. To have the right order of branching ratio Br{fi — )■ 67) as shown in the figure, we may have 
to consider, for example, the resonance effect. In this case two charged Higgs scalars are almost 
degenerate. For example, with the charged scalar mass ratio equals to 1.00006 ( green-solid line) 
we have the branching ratio i?r(/i — )■ 67) = 2.4 x 10~^^ for a value of A = 4.0 x 10~^°GeV~^ 
and the corresponding different /2 coupling constants are given as /J*^ = —4.55 x 10~^ and /J^ = 
(1.12 + 1 1.49) X 10-^ while, the / couplings are given by /"^ = 2.43 x 10"^ and = 6.20 x 10^^. 
We have seen the same order of magnitude appear once more for the same of charged scalar mass- 
ratio at a relative larger value of A. Here, we see for A = 7.05 x 10~^GeV~^, the branching ratio is 
Br{fi 67) = 1.60 X 10-12^ 




3e-09 6e-D9 9e-09 IJe-OS l.5e-0B 1.8e-0S 2.4e-(l 



A(GeV''f 




3^-09 6e-U9 9^-09 Lli-Oi 1.5e-0S l.Se-OS 2.1i^-0S 2.4e-0S 2.7e-C 
A IGtV"') 




3i;-09 6e-09 9e-09 1.2e-0S 1.5e-08 l.S^-OS 2.1^-OS 2.4e-0S 2.7e-C 
A |GeV''l 



FIG. 9: Variations of BR(/i ej) (left), BR(r ^7) (middle) and BR(r ej) (right) with the parameter 
"A" for different set of values (M^+/M^+). 



To check, should the same range of A and, hence, coupling constants will allow r-decay branching 
ratio to be within the experimental limit, we have shown the variation of both Br{T fi'-f) or 
Br{T — 67) on the middle and right panel of Fig.-9. Here, all the parameters are as mentioned above 
for the yU — )■ 67 discussion. We see that both the branching ratios can reach to the experimental 
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limit for the same range of A. For example, with the mass ratio equals to 1.00006 (green-solid line), 
we have the value of the branching ratio Br{T — )■ ^^7) = 1.73 x 10"^ and i?r(r — )■ 67) = 1.02 x 10~^ 
for a value of A = 7.05 x lO-^GeV"^ 



IV. FLAVOUR VIOLATING DECAYS INTO THREE CHARGED LEPTONS {k Ijkh) 



In this section we discuss the flavour changing charged trilepton decay modes. Recall that in our 
scenario, nonzero couplings are /'^'^, Z^"^, Z"^"^, /'^'^ from / matrix while from /2 matrix these are /J^ 
and 72'^ only. These nonzero couplings would induce /j — )■ Ijlkh decay with the exchange of neutral 
Higgs. 

At the tree level^, the matrix elements for the k — ?■ Ijlkk decay are given by 



Mh 



m 



H 



mi 



mi 



, cosa 
' t'sin/3 

sina 



rriE 



't'sin/? 

, 1 
' wsin/3 



mE + 



rriE 



cos(^^a) . 
V2smp 

V2smp 



. cosa cos(/3 — a) „ , ^ 

{—-^rriE ^ . ^ /2 Pl 

vsmp 

sina 



'wsin/3 



rriE + 



V2sm(3 
sin(/3 — a) 



f2)PL 



v^sin/3 



R 



' t'sin/3 



rriE 



y2sin/3 



f2)PL 



(26) 



The r decay width, thus, is given by 

r(r ^ //j) = 



3072 vr^ 



(27) 



where Dnj^ia are different contribution to the decay width from respective neutral Higgs scalars. 
Explicit expressions, for our model, one can read from Appendix-B. 

In case, two of the final fermions are identical, the antisymmetrization of the identical fermions 
and a symmetry factor 1/2 has to be considered in the formula which lead to an extra-factor 2 in 
the numerator. 



In Fig.-lO, variation of the branching ratio Brir — t- /i /i+yU ) and Brir — )■ e ) are shown 

on the upper row for a set of charged scalar mass ratio, m^+/m^+, ranging from 1.2 to 2.0, with 



set at 250 GeV. The corresponding variations of Brir — )■ /i 



e ' e 



and Brir e 



e ' e 



are shown in the lower row. On the top-left figure, we see that the branching ratio Brir^ — )• 
/i~/i"'"/i~) is approaching to the the present experimental limit < 1.8 x 10~^[6] for the variation 
of the parameter A upto 0.0002. The rest of the three graphs, within the same range of A, are 



^ For one loop contribution of trilepton decay modes of r-lepton see for e.g Mitsuda and Sasaki of Ref. [7]. We have 
also neglected very small contribution from the radiative r-decay processes, suprressed by a factor of olem-, into 
three charged leptons. 
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FIG. 10: Variations of the branching ratios of different r decay modes with the parameter "A" for different 
values of (iV4+/Af^+). 

showing the variation of the branching ratios for three other decay modes. For example, from the 
left-top figure we see that with the mass ratio equals to 1.4 ( red-dash line) we have the branching 
ratio Br{T~ — fi^jj^jj^) = 1.52 x 10"^ for a value of A = 1.65 x 10~"^GeV~^. So, in the same 
way one can read different coupling constants as before. With the same input parameter A, for the 
branching ratio Br{T~ — )■ e~e~^e~) the value is, from right bottom figure, 2.11 x 10~^^, as expected 
from the electron to muon mass-squared ratios. The corresponding values for Br{T~ — )■ e~/i+/i~) 
and Br{T~ — fi~e^e~) are given by 4.52 x 10~^° and 1.78 x 10~^^ respectively. In the analysis, we 
see that the Br{T /i/i/i) is the largest one, equals ~ 10^^ for a value of A ~ lO^^GeV"^. The 
same value of A and other input parameters we see that to other charged tri-leptons decay modes 
branching ratios are suppressed by the electron to muon mass ratio or it's higher power depending 
on the number of electron and muon in final states. 

The trilepton decay mode — > e"e"'~e"), at tree level, does not exist as it depends on the 
coupling constants /;^^. However, to add, we have to consider a scenario with a nonzero /l**^ if we 
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have any experimental evidence of muon decaying to three charged leptons process. 



V. CONCLUSIONS AND DISCUSSIONS 

Observation of neutrino oscillations, and hence neutrino masses, is one of the main evidence to 
look into beyond the standard model. Understanding of the smallness of neutrino masses requires 
new physics, such as loop induced masses as in the Zee model. These models do allow some pro- 
cesses which violate lepton flavour. There are strong experimental constraints on LFV interaction. 
Recently MEG collaboration has report a new upper bound of 2.4 x 10~^^ for — )■ 67 branching 
ratio. This improved upper bound may have implications on models for neutrino mass and mixing. 
On the other hand, the recent data from T2K, MINOS and Double Chooz provide some new infor- 
mation on the mixing angles in Vpmns that the last mixing angle 613 is non-zero at more than 2a 
level. The combined data analysis gives the confidence level at more than 3a. 

In the simplest Zee model, the mass matrix has all diagonal entries to be zero. This type of mass 
matrix has been shown to be ruled out. This is basically because that it cannot simultaneously 
have solution for |V^3| close to l/\/2 and \ Ve2\ close to as data require. 

In this article, it has been proposed that an interesting mass matrix can result if one imposes the 
requirement that no large hierarchies among the new couplings, that is, all nonzero f^^ and /g"' are 
of the same order of magnitude, respectively. This will lead to mass-matrix where all the elements 
are either proportional to the charged lepton mass or to it's square. A few important outcome in 
the neutrino section of the model are the foUowing:- 

• The model is constraint to have a non-zero ^13, which is in consistent with the present exper- 
imental data. 

• The best-fit value of our model predicts for 6*13 equals to 8.91° which is close to the recent 
data from Double Chooz. 

• Existence of solutions for non-zero CP violation with the Jarlskog parameter predicted in the 
range ±0.039, ±0.044 and ±0.048 respectively for a la, 2a and 3a ranges of neutrino masses 
and mixing angles. However, lifting the constraints on ^13 the above respective ranges become 
±0.049, ±0.053 and ±0.056. 

For the inverted hierarchy only, the best-fit values of the mass matrix parameters are: 



and, the corresponding output for the mixing angles and mass-squared differences are given as 



We have also discussed if these best-fit values satisfy different constraints for the lepton flavour 
violating processes. In that direction we have considered the radiative decay processes like — )■ 67 



X = 0.255, y = 4.100, z = 1.790, a = 0.017 eV, 6 = 180°, 



follows 



sin^^ia = 0.3163, sin^^aa = ( 
W = 7.51 X lO-^eV^ Am^ 



0.4033, sin2^i3 = 0.0256, 
^ = -2.36 X 10~Wl 
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and r /u(e)7 processes and the decay modes to three charged leptons r~ — )■ fi~fi~^fi~, t~ — )■ 
e~/i~^yU^, — 7- fi"e~^e~ and r~ — j- e~e"'"e~ processes. Here, we look into that range of parameter 
space which will lead to the experimental limit of these decay modes, for example, the constraints 
set by the MEG collaboration for the /i — )■ 67 process. 

To have a better understanding on how different parameters contribute to the branching ratios 
we have discussed different cases - when either of the coupling constants /2 or / has a dominant 
contribution and when both the couplings / and /2 have comparable contributions to different n 
and r-radiative decay modes. 

We have extensively shown how these different cases correspond to various ranges of the param- 
eter A, and hence the set of non-zero coupling constants, to have a branching ratio close or within 
the experimental limit. 

There are a pair of charged Higgs bosons, two CP-even and one CP-odd neutral Higgs bosons 
in the Zee model. The Higgs sector is just enriched by one more charged scalar component than 
the Two Higgs doublet model. The collider phenomenology of Higgs sector in the Zee model is 
studied in various literature, for example [20]. A very well known way to analyse this model at 
LHC, is the pair production of charged and neutral Higgs via Drell-Yan process. Different decay 
modes of these Higgs bosons will lead to multi-lepton plus missing energy topology in the final 
states. However, depending on the coupling constants and mass of the neutral Higgs boson various 
channel, like decaying to di-photon, will be modified. How different channels are affected due to 
our model specific coupling constant parameter space will be separately discussed in future. 

Within the allowed non-zero parameter space, the only possible lepton flavour violating decay 
modes in our model are /x — )■ 67, r — )■ /i(e)7, r — )■ ee) and r — e(/i/i, ee). Collider search of 

lepton-fiavour violating decay modes is a clear signal to hunt for beyond the SM. A large number 
of past as well as recent articles have already studied this LFV. In our analysis, we see that the 
Br(T — /i/i/i) is the largest one, equals ~ 10^^ for a value of A ~ lO^^GeV^^. Not only a relatively 
large branching ratio but muons are favourable for their clear signal at LHC. We would thus be 
interested here for the decay channel r — ?■ fififi mediated via different neutral scalars. At LHC, the 
main contribution to the production of r-lepton via the heavy meson B or D decays and weak W 
and Z gauge bosons. Recently, LHC acquired a total integrated luminosity of 5 fb~^ per year. We 
will consider this analysis in a future work [21] in detail. 
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Appendix A: Higgs Scalars Contribution to Radiative LFV Decays 



In this appendix, we write down the exphcit expression of Cl and Cr due to the charged scalars, 
hi,h2, and neutral scalars, H, h, a. These two factors Cl and Cr can be written as 



Cl = Y^{Q5(c,+ + C7,+)-g;(c^ + c, + a)} 



(Al) 



where, Qs = —1 is the charge of /ij and Q/ = — 1 is the charge of charged lepton, while Cr = C\^. 
The matrix is given by 



C 



H 



Ca 



cosa ^ cos(/3-a) ft \ p \ I cosa_ cos(p-a) ,■ \ 

J. m (^°^m„ cos(/3-Q!) r\ p / \ / cosa cos(/3-a) ff 



cos(/3— a) 



I / cosa_ cos(ff-a) f \ T? ( -y \ I cosa ^ cos(,p-Qj j. 



cos(/3— a) 



2 ' 
J? 



sin(/3-a) ft \ p \ / sing ^ , sm(p-a) r \ 



sin(/3— a) 



_L ^ / sing _|_ sin(^-g) f \ rp ( ^ \ I sing , smyp-g; 



sin(/3-g) 



I / sing I sin(/3-g) f \ p / \ / sing , sin(/3-g) f 

+ [^^E + -y^h) mE^3{Zh) [^ruE + ^/^h 

+ ^E [^ruE - 75^-^/2) F,{za) [^rriE - 
+ {■^'^E - ^/2) msF3(^,) (^ms - 



2 ' 



m 



2 ■ 



(A4) 



(A5) 



(A6) 



In the above equations Fa{zi) is a diagonal matrix Fa{zi) = diag{Fa{zi), Fa{z2), Fai^z^)), where 
Zi = ml/m^g with rrti the fermion "i" mass and scalar "S" in the loop, respectively. 

The functions Fi{z)^F2{z) and -^3(2;) are given by 



F,{z) 



z"^ — 5z — 2 zlnz 



12(2-1)3 2(2-1)4 ' 

22;2 -)- 5^; — 1 z^hiz 

12(2-1)3 ~ 2(Z- 1)4 
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To be noted here, in the hmit z both the functions Fi{z) and F2{z) do not any more depend 
on the argument of the functions, and take finite values —1/6 and —1/12 respectively. On the 
other hand, the characteristic of the function F^lz) is different and goes as In 2; in the limit ^ — ?■ 0. 
However, the final expressions for Cl,r are finite in the limit z ^ 0. 



Appendix B: Neutral Higgs scalar contribution to r — )• hhh 



Here, we will explicitly show different neutral scalar contribution to the lepton flavour violating 
trilepton decay mode of r-decays in our model. With the non-zero / and /2 we are considering, we 
find only, r — )■ ee) and r — )■ e(yU/2, ee). 

The matrix elements are of the form: M = DalljPnT with given by for different processes 
and different Higgs a contributions 

^ 1 cosa cos(/3 - a) 

T-^fififi: Dh = 2 ^ — /K • o ) 

mjf vsmp V2sm/3 

1 sing sin(/3 - a) 

ml vsmfj V2sm/3 



ml vsin/S ^/2sm(3 
1 cosa cos(/3 — a) 



re* 



mjj vsm(3 V2sm/3 
_ 1 sina sin(/3-a) 
wsm/^ V2sm/3 

n — ^ I I rn f^^* 

wsm/i( V2sm/3 

1 cosa cos(/3 - a) 

r ^ /xee : = 2 ^ — /k • ^ "^^/a , 

i;sm/^ V2sm/3 

1 sina sin(/3 - a) 

t;sm/3 V2sm/3 

mi vsmp y/2sm(3 
_ ^ 1 cosa cos(/3 - a) 

mjj vsm.[i V2sm/3 
_ 1 sina sin(/3-a) 
f sm/? V2sm/3 

1 1 1 

n — rn r''* 

mi vsm.p V2sm/3 
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